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ABSTRACT 


This  thesis  focuses  on  the  modeling  and  simulation  of  the  performance  of  M-ary 
modulation  techniques.  In  addition  to  the  most  popular  used  coherent  modulation 
schemes,  forward  error  correction  (FEC)  coded  signals,  and  fading  channels  are  also 
considered.  The  channel  is  modeled  for  the  two  cases  of  fading  and  no  fading,  both  with 
additive  white  Gaussian  noise  (AWGN).  The  effect  of  barrage  noise  interference  is  also 
considered.  Binary  phase-shift  keying  (BPSK),  16-PSK,  and  16-QAM  (quadrature 
amplitude  modulation)  communication  systems  are  simulated  and  analyzed.  For  BPSK, 
both  hard  decision  decoding  (HDD)  and  soft  decision  decoding  (SDD)  are  considered. 
For  16-PSK  and  16-QAM,  only  HDD  is  examined.  Simulation  results  and  analytical 
results  are  compared.  The  results  show  that  the  modeling  and  simulation  in  SystemView, 
are  in  excellent  agreement  with  the  analytical  results. 
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EXECUTIVE  SUMMARY 


The  purpose  of  this  thesis  is  to  model  and  simulate  the  most  commonly  used 
coherent  systems,  such  as  binary  phase-shift  keying  (BPSK),  16-PSK,  and  16-QAM 
(quadrature  amplitude  modulation)  for  different  kinds  of  noisy  and  fading  channels  in 
order  to  develop  a  mechanism  to  evaluate  the  performance  of  these  systems  for 
conditions  that  do  not  easily  yield  to  analysis.  This  is  important  because  these  coherent 
modulations  are  used  in  many  communication  systems  including  the  sub-carriers  in 
Orthogonal  Frequency  Division  Multiplexed  communications  systems  such  as  IEEE 
802.11a,  g,  and  n  WLANs  and  the  Joint  Tactical  Radio  System  Wideband  Networking 
Waveform.  Therefore,  in  order  to  predict  the  performance  of  these  systems,  simulations 
typified  by  the  simulations  in  this  thesis  are  necessary.  SystemView  provides  a  method  to 
simulate  the  performance  of  communication  systems  that  are  otherwise  difficult  to 
analyze.  In  this  thesis,  coherent  modulation  schemes  with  forward  error  correction  (FEC) 
coded  signals  transmitted  over  fading  channels  are  simulated.  The  channel  is  modeled  for 
both  fading  and  no  fading  with  additive  white  Gaussian  noise  (AWGN).  The  effects  of 
barrage  noise  interference  are  also  considered.  BPSK,  16-PSK,  and  16-QAM 
communication  systems  are  simulated  and  analyzed.  For  BPSK,  both  hard  decision 
decoding  (HDD)  and  soft  decision  decoding  (SDD)  are  considered.  For  16-PSK  and  16- 
QAM,  only  HDD  is  examined. 

In  order  to  verify  the  modeling  and  simulation  procedure,  the  performance  of 
coherent  systems  with  no  coding  was  examined  first.  Next,  the  performance  of  coherent 
systems  with  error  correction  coding  was  examined.  Both  hard  and  soft  decision 
decoding  was  considered  for  binary  systems,  but  only  hard  decision  decoding  was 
considered  for  non-binary  systems.  Finally,  fading  channels  were  also  considered. 

Two  major  problems  were  found  during  the  modeling  and  simulation.  First,  the 
total  number  of  samples  is  restricted  to  231  (=  2,147,483,647)  for  all  simulation  loops, 
which  reduces  accuracy  for  high  signal-noise-ratios  (SNRs).  Next,  too  much  CPU 
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resource  is  required  to  conduct  each  simulation.  These  restrictions  are  an  artifact  of  using 
the  student  version  of  SystemView  and  presumably  would  not  be  as  much  of  a  drawback 
with  the  commercial  version. 

The  performance  of  the  coherent  systems  examined,  both  with  and  without  coding, 
all  agreed  very  well  with  analytical  results.  Since  the  simulations  developed  have  been 
verified  for  the  AWGN  case,  they  can  now  be  modified  to  examine  cases  that  are  not 
possible  to  evaluate  by  analysis. 


I. 


INTRODUCTION 


A.  THESIS  OBJECTIVE 

Wireless  local  area  networks  (WLANs)  are  increasingly  used  for  both  commercial 
and  military  applications.  The  objective  of  this  thesis  is  to  investigate  the  performance  of 
systems  that  use  different  combinations  of  modulation  order  and  code  rates  of  forward 
error  correcting  codes  such  as  those  typically  used  in  WLANs  implemented  with 
orthogonal  frequency-division  multiplexing  (OFDM).  This  thesis  focuses  on  modeling 
and  simulation  of  the  performance  of  M- ary  modulation  with  forward  error  correction 
(FEC)  coded  signals  for  both  channel  fading  and  no  channel  fading  with  additive  white 
Gaussian  noise  (AWGN).  Both  hard  decision  decoding  (HDD)  and  soft  decision 
decoding  (SDD)  are  considered.  Since  the  analysis  with  SDD  is  not  exact  when  using 
non-binary  modulation  and  a  binary  code,  it  is  important  to  verify  approximate  analytical 
approaches  by  simulation. 


B.  RELATED  RESEARCH 

Wireless  communications  have  ushered  in  a  new  era.  Numerous  studies  have 
investigated  the  performance  of  OFDM  for  many  different  fading  channels,  such  as 
OFDM  over  frequency-nonselective,  fast  Rayleigh  and  Ricean  fading  channels  in  [1]  and 
[2],  frequency-selective,  Rayleigh  fading  channels  in  [3],  frequency-selective,  slow 
Ricean  fading  channels  in  [4],  frequency-selective,  slow  Nakagami  fading  channels  in 
[5],  and  frequency-selective,  slow  Nakagami  fading  channel  with  pulsed  jamming  in  [6]. 
Theses  papers  derived  analytical  expressions  for  the  different  channel  characteristics. 
These  are  just  a  few  of  many  related  studies.  Distinct  from  the  above  works,  in  this  thesis, 
we  focus  on  the  modeling  and  simulation  of  the  performance  of  M- ary  modulation  with 
FEC  coded  signals  over  fading  channels.  The  channel  is  modeled  for  the  two  cases  of 
fading  and  no  fading,  both  with  AWGN.  Binary  phase-shift  keying  (BPSK),  16-PSK,  and 
16-QAM  (quadrature  amplitude  modulation)  communication  systems  are  simulated  and 
analyzed. 
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C.  THESIS  ORGANIZATION 

After  the  introduction,  the  thesis  is  organized  into  five  remaining  chapters. 
Following  an  introduction,  Chapter  II  introduces  the  fundamental  background  of  the 
digital  communication  system,  including  the  modulation  techniques.  Fading  channels  are 
discussed  in  Chapter  III.  Forward  error  correction  (FEC)  is  also  introduced  in  this  chapter. 
Next,  Chapter  IV  explores  the  differences  between  simulation  results  and  the  analytical 
results.  The  use  of  SystemView  software  to  simulate  the  system  under  investigation  is 
also  discussed  in  this  chapter.  Finally,  Chapter  V  is  a  summary  and  recommendation  for 
future  work. 


D.  REQUIRED  SOFTWARE 

SystemView  by  Elanix,  version  5.0  was  used  to  build  and  conduct  the  simulations. 
MATLAB  version  7.1.0  (R 14)  was  used  to  generate  the  diagrams. 
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II.  DIGITAL  COMMUNICATION  SYSTEM 


Digital  Communication  has  entered  our  daily  lives  in  many  different  ways,  for 
example,  the  radio  and  television  sets  in  our  living  rooms  and  the  networks  around  us.  In 
this  chapter,  the  fundamental  background  of  digital  communication,  including  various 
modulation  techniques  and  fading  channels  are  described. 


A.  OVERVIEW 

Digital  communication  can  be  simply  thought  of  as  a  series  of  processes  from 
information  input,  encoding,  modulating,  demodulating,  decoding,  and  information 
recovery.  The  basic  digital  communication  model  is  described  in  Figure  1 . 

With  digital  binary  communication  systems,  one  of  two  unique  signals,  say  (t) 
or  s2(t) ,  which  represent  a  bit  “1”  or  a  bit  “0”,  are  transmitted  each  bit  duration. 

Normally,  the  signals  may  be  represented  as  a  time  sequences  of  symbols,  where  there 
are  M  distinct  symbols,  and  each  symbol  represents  q  bits  of  information.  Consequently, 
the  relation  between  symbol  alphabet  size  M ,  and  q  bits  of  information  can  be  expressed 
as, 

q  =  l°g2  M  (2.1) 


or 


2 q  =M 


(2.2) 
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Figure  1.  Digital  Communication  Model  (After  Ref.  7) 


Shannon’s  theorem  says  that  for  an  arbitrarily  small  probability  of  error,  the 
maximum  possible  bandwidth  efficiency  is  limited  by  the  noise  in  the  channel  [8],  For  a 
channel  with  AWGN,  the  capacity  can  be  stated  as  [9] 


C  =  W  log2 


n  J 


(2.3) 


or 


(2.4) 


C  =  W  log2 


^  E  C  ^ 
1  +  -  fcC 


N0Wj 


where  C  is  the  channel  capacity  in  bits/s,  Pc  is  the  average  received  signal  power,  Pn  is 
the  average  received  noise  power,  Eb  is  the  average  received  signal  energy  per  bit,  W  is 
the  noise  equivalent  system  bandwidth  in  Hz,  and  N0  /  2  is  the  noise  power  spectral 
density  (PSD)  in  W/Hz. 
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For  a  certain  Ebl  N0,  Shannon’s  theorem  states  that  if  C>Rb,  where  Rb  is  the  bit 
rate,  then  it  is  possible  to  achieve  error  free  communication.  On  the  contrary,  if  Rb  >  C , 

then  error  free  communication  is  not  possible.  For  a  general  expression,  Equation  (2.4) 
was  modified  as  follows: 


Ejl 

N0 


(2C/W  -1) 


C/W 


(2.5) 


This  equation  is  typically  plotted  on  the  bandwidth  efficiency  plane  ( Rb/W 
versus  Eb/  N0)  [11]  as  shown  in  Figure  2.  If  Rb  is  the  data  rate  in  bits  per  second  and  W 
is  the  bandwidth  occupied  by  the  modulated  radio  frequency  (RF)  signal,  then  bandwidth 
efficiency  ijH  is  expressed  as  [10] 


(2.6) 


According  to  this  equation  above,  there  are  two  different  regions  on  this 
bandwidth-efficiency  plane.  The  bandwidth-limited  region  has  rjB>  1 .  The  power-limited 

region  has  rjB<  1 . 
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R/W  (bit/s/Hz) 


B.  MODULATION  TECHNIQUES 

1.  Binary  Phase-Shift  Keying  (BPSK) 

In  binary  phase-shift  keying  modulation,  it  represents  a  bit  “1”  or  a  bit  “0”  by 
changing  the  phase  of  the  waveform,  either  zero  or  n  radians  (180°),  respectively.  The 
waveform  with  this  sudden  phase  change  is  shown  in  Figure  3. 
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A  general  transmission  signal  is 

s(t)  =  V2 Ac  cos  [2 n fct  +  Gd  (t)  +  60  ]  (2.7) 

where  2n fc  is  the  carrier  frequency  in  radians/s,  6d (t)  represents  the  data,  and  0o 
is  the  fixed  reference  phase.  Frequently,  0d{t)  =  0  is  used  to  represent  bit  “1”,  while 
6d  (t)  =  7i  is  used  to  represent  bit  “0”. 


Baseband  Data 


1 


0  0 
_ I _ 


1 


Binary  PSK 
waveform 


Figure  3.  Data  and  waveform  format  for  BPSK  (After  Ref.  12) 


2.  Multi-level  (M-ary)  Phase  and  Amplitude 
a.  M-ary  Phase  Shift  Keying  (MPSK) 

On  the  bandwidth  limited  region,  a  frequently  used  modulation  format  is 
M- ary  phase-shift  keying  (MPSK).  The  transmission  signal  of  the  ith  symbol  is 


sit)  =  V2 Ac  cos 


2*/c,+M!tz2)+6, 

M 


(2.8) 


where  m  =  1,  2,  ...,  M. 

Unlike  BPSK  which  only  utilizes  the  in-phase  component,  MPSK  uses  both  in- 
phase  and  the  quadrature  components. 

b.  M-ary  Quadrature  Amplitude  Modulation  (MQAM) 

The  constellations  for  both  16QAM  and  16PSK  are  shown  in  Figure  4.  In 
MQAM,  amplitude  and  phase-shift  keying  are  combined  to  transmit  several  bits  per 
symbol  (in  this  case  q  =  4  bits/symbol). 
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16QAM 


16PSK 


Figure  4.  The  constellation  of  16QAM  and  16PSK  (After  Ref.  12) 


In  this  chapter,  digital  communication  systems  were  introduced  and  modulation 
schemes  such  as  BPSK,  MPSK  and  MQAM  were  discussed.  In  the  next  chapter,  the 
primary  sources  of  channel  interference  are  considered  and  error  correction  coding  is 
introduced. 
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III.  FADING  CHANNELS  AND  ERROR  CORRECTION  CODING 


Additive  white  Gaussian  noise  (AWGN)  is  always  present  in  all  communication 
systems  and  affects  each  transmitted  symbol  independently  [11].  In  addition  to  AWGN, 
the  focus  of  this  chapter  is  on  fading  channels  and  error  correction  coding. 


A.  FADING  CHANNELS 

Reflection,  diffraction,  and  scattering  are  the  main  physical  factors  that  affect  the 
signals  to  be  received.  Due  to  these  factors  radio  frequency  channels  are  extremely 
random  and  are  not  easily  analyzed.  Many  wireless  communication  channels  do  not  have 
a  line-of-sight  (LOS)  signal  path.  When  there  is  no  LOS,  the  signal  is  only  received  by 
multipath ;  that  is,  there  are  multiple  signal  paths  from  the  transmitter  to  the  receiver  as  a 
result  of  reflection,  diffraction,  and  scattering  of  the  original  signal.[ll]  Multipath  is 
illustrated  in  Figure  5. 


Figure  5.  Multipath  propagation 


Multipath  can  cause  large  fluctuations  in  the  received  signal’s  amplitude,  phase, 

and  angle  of  arrival.  Figure  6  represents  an  overview  of  how  fading-channel  can  be 

categorized.  Large-scale  fading  relates  received  signal  strength  for  large  changes  in 

transmitter  and  receiver  separation  distance  (usually  much  bigger  than  signal’s 

wavelength).  Small-scale  fading  (changes  in  separation  distances  on  the  order  of  a  half 

wavelength)  is  used  to  describe  the  rapid  fluctuation  of  the  amplitude  of  the  received 

signal  over  a  short  period  of  time,  and  large-scale  fading  effects  may  be  ignored.  [11] 
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Figure  6.  Fading-channel  manifestations  (From  Ref.  1 1) 


B.  SMALL-SCALE  FADING 

The  distance  between  small-scale  fades  is  on  the  order  of  A  /  2 ,  as  shown  in 
Figure  7.  As  Figure  6  shows,  the  two  main  mechanisms  that  cause  small-scale  fading  are 
time  spreading  of  the  signal  due  to  multipath  and  time  variance  of  the  channel  due  to 
Doppler  spread. 


Power/ 

dBm 


Small  Scale 
Fading 

Large  Scale 
Fading 


Distance 


Figure  7.  Large  and  small  scale  fading  (From  Ref.  13) 


10 


1.  Time  Spreading  of  the  Signal  Due  to  Multipath 

As  shown  in  Figure  6,  time  spreading  of  the  signal  will  can  be  due  to  either 
frequency-selective  fading  or  frequency  nonselective  fading  (flat  fading).  The  coherence 
bandwidth  f0  and  the  maximum  excess  delay  Tm  are  two  factors  used  to  describe  fading. 

The  reciprocal  relationship  between  f0  and  Tm  is  given  by  [1 1] 

(3-1) 

1  m 

where  the  time  Tm  represents  the  maximum  excess  delay,  or  the  difference  of  arrival 
times  of  the  first  and  the  last  received  signal  components. 

Typically,  the  received  signal  consists  of  several  different  multipath  components,  which 
causes  the  average  received  power  S(r)  to  present  multiple  isolated  peaks,  as  illustrated 
in  Figure  8.  [1 1] 

Six) 


Maximum  excess  delay 

Figure  8.  Multipath  intensity  profile  (After  Ref.  1 1) 

The  multipath  intensity  profile  can  change  significantly  for  channels  with  the  same  value 
of  Tm .  A  general  parameter,  the  root-mean-square  (rms)  delay  spread,  was  used  instead  of 

Tm ,  defined  [11] 
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/T  “2 

(7t=\It  -t 


(3.2) 


where  r  is  the  mean  excess  delay  and  r2  is  the  second  moment  of  r  . 

There  is  no  fixed  relationship  between  /0 ,  the  coherence  bandwidth,  and  crr ,  the 

root-mean-square  (rms)  delay  spread.  If  the  coherence  bandwidth  is  defined  as  the 
frequency  interval  over  which  the  frequency  correlation  function  has  a  correlation  of  at 
least  0.9,  the  coherence  bandwidth  is  approximately  [10] 


/„ 


1 

50crr 


(3.3) 


If  the  frequency  correlation  function  is  above  0.5,  then  the  coherence  bandwidth  is 
empirically  obtained.  [14] 


(3.4) 


a.  Frequency-Selective  Fading 

For  the  frequency-selective  case,  if  the  coherence  bandwidth  f0  is  less 
than  the  symbol  rate  1  /  Ts ;  that  is  [  1 1  ] 

(3-5) 


which  means 


T,  >  T„  (3.6) 

where  the  symbol  rate  l/Ts  is  equal  to  the  signaling  rate.  Figure  9(a)  shows  that  various 
spectral  components  of  the  transmitted  signal  will  be  affected  differently. 


b.  Flat  Fading 

Frequency-nonselective  or  flat-fading  degradation  appears  whenever  [11] 
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A»r 


(3.7) 


or,  equivalently, 


T  <TC 


(3.8) 


Figure  9(b)  illustrates  that  all  of  the  spectral  components  of  the  transmitted  signal  are 
affected  in  approximately  the  same  way. 


Spectral 

density 


Transmitted 


(a)  Typical  frequency-selective  fading  case  (fo  <  W ) 


Spectral 

density 


Channel 


Figure  9.  Typical  cases  for  frequency-selective  and  flat-fading  (From  Ref.  11) 

2.  Time  Variance  of  the  Channel  Due  to  Doppler  Spread 

There  are  two  parameters  used  to  describe  the  time- variant  nature  of  the  channel: 
Doppler  spread  fd  and  coherence  time  T0 . 


The  magnitude  of  the  frequency  shift  fd  is  given  by  [1 1] 


(3.9) 
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where  V  is  relative  velocity  and  X  is  the  signal  wavelength.  When  the  transmitter  and 
receiver  comes  to  each  other,  fd  is  positive,  and  when  they  move  away  from  each  other, 

fd  is  negative. 


The  Doppler  spread  fd  and  the  coherence  time  7)  are  reciprocally  related.  [11] 


(3.10) 


a.  Fast  Fading 

A  channel  is  said  to  be  fast  fading  if  the  symbol  rate  1/7).  (approximately 
equal  to  the  bandwidth  IF)  is  less  than  the  fading  rate  1/7)  (approximately  equal  to  fd); 
that  is,  fast  fading  is  characterized  by  [1 1] 

W<fd  (3.11) 

or 

TS>T0  (3.12) 

In  practice,  fast  fading  only  occurs  for  low  data  rates  [11], 


b.  Slow  Fading 

In  the  frequency  domain,  this  implies  that  the  Doppler  spread  of  the 
channel  is  much  less  than  the  bandwidth  of  the  baseband  signals.  That  is  [1 1] 

W>fd  (3.13) 


or 


(3.14) 
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3.  Summary  of  Small-Scale  Fading 

Depending  on  the  relation  between  the  signal  parameters  (such  as  bandwidth, 
symbol  period,  etc.)  and  the  channel  parameters  (such  as  rms  delay  spread  and  Doppler 
spread),  different  transmitted  signals  will  experience  different  types  of  fading.  Signal 
time  spreading  (signal  dispersion)  and  the  time-variant  nature  of  the  channel  are  two 
mechanisms  that  result  in  small-scale  fading.  Figure  10  summarizes  these  two  small-scale 
fading  mechanisms.  [11] 
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Figure  10.  Summary  of  Small-scale  fading  (After  Ref.  11) 


C.  RAYLEIGH  AND  RICEAN  CHANNELS 
1.  Rayleigh  Fading  Channel 

In  the  Rayleigh  fading  case,  there  is  no  LOS  signal  transmitted  between 
transmitter  and  receiver  and  all  of  the  received  signal  power  is  due  to  multipath.  For 
Rayleigh  fading  channels  (a  special  case  of  the  Ricean  fading  channel),  the  received 
amplitude  is  modeled  as  a  Rayleigh  random  variable  with  the  probability  density  function 
(pdf)  given  by  [14] 


p(r)  =  —Qxp 


2  \ 


2<t 


,  (0  <  r  <  oo ) 


(3.15) 
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where  a  is  the  rms  value  of  the  received  voltage  signal  before  envelope  detection  and 
cr2  is  the  time-average  power  of  the  received  signal  before  envelope  detection. 


2.  Ricean  Fading  Channel 

When  there  is  a  dominant  signal  component  present,  such  as  a  LOS  signal,  the 
effect  of  a  LOS  signal  arriving  with  many  weaker  multipath  signals  gives  rise  to  the 
Ricean  distribution.  The  pdf  of  the  Ricean  distribution  is  given  by  [14] 


p(r)  = 


(r2  +  A2) 

L 

(Ar) 

1 

<N 

b 

<N 

_ l 

0 

A 

(A  >  0,  r  >  0) 


(3.16) 


The  parameter  A  denotes  the  peak  amplitude  of  the  dominant  signal  and  70(»)  is  the 
modified  Bessel  function  of  the  first  kind  and  zero-order. 

The  Ricean  distribution  is  often  described  in  terms  of  a  parameter  C,  (Ricean  factor)  [14], 
which  is  defined  as  the  ratio  between  the  deterministic  LOS  signal  power  and  the  non- 
LOS  multipath  signal  power,  given  by  £  =  A1  /(2cr2) ,  or  in  dB , 

A2 

C(dB)  =  10  log — -dB  (3.17) 

2(7 

As  A  — >  0,  C,  — >  -oo (dB) ,  and  as  the  dominant  path  decreases  in  amplitude,  the 
Ricean  distribution  degenerates  to  a  Rayleigh  distribution. 


D.  SUMMARY  OF  MULTIPATH  FADING  CHANNELS 

The  characteristic  of  the  channel  was  described  and  multipath  fading  channels 
were  discussed.  The  mechanisms  and  types  of  small  scale  fading  channels,  frequency- 
selective  fading,  flat  fading,  fast  fading,  and  slow  fading  were  also  introduced.  In  this 
thesis,  two  widely  used  models  for  fading  channels,  the  Rayleigh  fading  channel  and  the 
Ricean  fading  channel,  are  assumed.  In  the  next  section,  error  control  coding  is  discussed. 
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E.  ERROR  CORRECTION  CODING 

Error  correction  coding  has  made  digital  communications  practical  and  popular. 
By  using  the  error  correction  coding  method,  digital  systems  can  easily  detect  and  correct 
errors.  All  communication  systems  should  be  designed  to  combat  the  degenerative  effects 
of  the  communication  channel  such  as  noise,  interference,  and  fading. 

There  are  two  basic  error  control  strategies:  automatic  repeat  request  (ARQ)  and 
forward  error  correction  (FEC)  coding.  Forward  error  correction  coding  consists  of 
adding  a  certain  number  of  redundant  bits  to  the  actual  data  bits  in  a  particular  pattern 
such  that  recovery  of  the  actual  data  bits  is  enhanced.  [1 1]  In  this  thesis,  FEC  is  used. 


1.  Forward  Error  Correcting  (FEC)  Coding 

Forward  error  correction  (FEC)  coding  consists  of  adding  a  certain  number  of 
redundancy  bits  to  the  actual  data  bits  in  a  specified  pattern  such  that  recovery  of  the 
actual  data  bits  is  possible.  The  two  primary  types  of  forward  error  correction  codes  are 
block  codes  and  convolutional  codes. [1 1]  In  this  thesis,  convolutional  codes  are  used. 

In  a  system  utilizing  FEC  coding,  for  every  k  information  data  bits,  n  coded  bits 
are  transmitted  where  n  >  k ,  the  relation  is, 

nTbc=kTb  (3.18) 

where  Th<  is  the  duration  of  a  coded,  or  channel,  bit  and  Tb  is  the  duration  of  a  data  bit. 
Therefore, 

Tlc  =-Tb=rTb  (3.19) 

c  n 

Kc=—  (3-20) 

r 


where  the  coded,  or  channel,  bit  rate  is  Rh  =  1  /  f  and  r  =  k/  n  is  the  code  rate. 

Since  r<  1,  the  coded  bit  rate  is  higher  than  the  uncoded  bit  rate;  and  for  a  constant  data 
bit  rate,  the  addition  of  FEC  coding  increasing  the  bandwidth  by  the  factor  n/k  =  l/r . 
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[11]  Also,  the  average  energy  per  coded  bit  is  less  than  the  average  energy  per  uncoded 
bit;  and  for  a  fixed  average  energy  per  data  bit,  the  addition  of  FEC  coding  increases  the 
probability  of  channel  bit  error.  [11] 

The  average  transmitted  power  is  the  same  whether  coded  or  uncoded  bits  are 
transmitted: 

EbRbc=EbRb  (3.21) 

Therefore, 

Ebc=rEb  (3.22) 


Coding  gain  is  defined  as  the  difference  in  the  signal-to-noise  ratios  required  by  a 
FEC  coded  communication  system  and  the  same  uncoded  communication  system  to 
achieve  a  specific  probability  of  bit  error: 


GdB  = 


^  e  ^ 

V  J  uncoded, r. 


f  E  ^ 

n b 

.  j\f  . 

^  ^  x  coded  j  i 


(3.23) 


2.  Hard  Decision  Decoding  (HDD)  and  Soft  Decision  Decoding  (SDD) 

In  the  process  of  decoding,  two  different  types  of  “decisions”  can  be  made  in 
order  to  recover  the  original  information.  A  “hard  decision”  can  be  made  in  which  a 
received  bit  is  demodulated  as  being  either  a  binary  one  or  zero.  This  determination  is 
then  sent  to  the  Viterbi  decoding  algorithm  for  decoding  of  the  actual  information  bits. 

A  “soft  decision”  may  also  be  made  in  which  the  decoder  allows  for  varying 
degrees  of  certainty  that  a  received  bit  is  a  one  or  a  zero.  This  quantization  enables  the 
decoder  to  classify  a  received  bit  as  a  “strong”  binary  one  or  zero  or  marginally  declaring 
the  bit  a  binary  one  or  zero.  This  approach  improves  the  ability  of  the  decoder  to 
accurately  estimate  a  received  bit’s  identity.  [11] 

In  this  chapter,  fading  channels  and  error  correction  coding  were  both  discussed. 
In  the  next  chapter,  the  differences  between  simulation  results  and  the  analytical  results 
are  compared. 


18 


IV.  COHERENT  SYSTEM  SIMULATION  IN  SYSTEMVIEW 


A.  PERFORMANCE  OF  COHERENT  SYSTEMS  WITH  NO  CODING 

1.  Coherent  Detection  (BPSK,  MPSK,  A/QAM) 

In  this  subsection,  the  probability  of  bit  error  for  binary  phase-shift  keying 
(BPSK),  M- ary  phase-shift  keying  (MPSK),  and  M- ary  quadrature  amplitude  modulation 
(MQ AM)  waveforms  transmitted  over  slow,  flat  fading  channels  is  listed.  The  following 
equations  and  tables  are  summarized  after  [15]. 

For  each  of  these  signaling  techniques,  the  conditional  probability  of  bit  error  is  of 
the  form  [15] 

PM  =  ^Q(^)  (4.i) 

where  Q{ •)  is  the  Q-function  and  a  ,  b  ,  and  q  are  constants  that  depend  on  the 
modulation  type.  The  constants  a ,  b ,  and  q  for  BPSK,  MPSK,  and  MQAM  are  listed  in 
Table  1.  For  MQAM  with  q  even,  a  square  constellation  is  assumed,  while  for  MQAM 
with  q  odd,  a  rectangular  constellation  is  assumed.  For  example,  in  Figure  11,  M  =  8 
(q= 3,  blue  line)  yields  a  rectangular  constellation  and  M  =  16  (q= 4,  red  line)  has  a  square 
constellation. 
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Table  1 .  Modulation  dependent  constant  for  Equation  (4.1)  (After  Ref.  15). 


Modulation 

q 

a 

b 

BPSK 

l 

i 

2 

QPSK 

2 

2 

1 

MPSK 

log  2M 

2 

2 -sin  2{n !  M') 

MQAM 

4,6,8,- 

4-(l-2^/2) 

3  /(2q  - 1) 

MQAM 

3,5,7,- 

4 

3  !(2q  - 1) 

2.  Model  the  Coherent  System  (BPSK,  MPSK,  MQAM)  in  SystemView 
a.  Simulation  of  BPSK  in  System  View 

Baseband  BPSK  is  the  best  binary  modulation  technique  in  terms  of  BER 
for  channels  with  AWGN.  In  SystemView,  a  PN  sequence  generator  is  used  to  simulate 
the  data  source.  Figures  12  and  13  show  the  simulation  results  for  baseband  BPSK,  which 
shows  good  agreement  with  the  analytical  result  in  equation  (4.1). 
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Bit  error  probability,  P, 


s 

■ 


Performance  of  BPSK  in  AWGN 
with  no  channel  fading. 


tfc — iflM 


Final  Value  Sink  10: 

Loop  N  Time  Value 

1  19080  19  0800Q00000000e-3  7.8G1G3522012579e-2 

2  57735  57.7350000000000e-3  3. 880481 179GGG28e-2 

3  1 78432  1 78. 432000000000e-3  1 . 242781 51 072520e-2 

4  80G32G  80G  326000000000c  3  2.388938260279G0e  3 

5  87591 38  8. 7591 3800000000  1 . 88G1 2531 01 41 87e-4 

G  758G8001  75.8680010000000  3.G209822240G1 82e-G 
7  14297G8G4  142.97G8G4500000  1.49011G14158923e-8 


Figure  12.  The  simulation  of  baseband  BPSK  in  SystemView. 


Performance  of  simulated  baseband  BPSK  in  AWGN  without  coding 


10 1 
10'2 
10'3 
10'4 
10'5 
10* 
10'7 
10* 


Figure  13.  The  differences  of  performance  between  simulated  and  theoretical  BPSK. 
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b.  Simulation  of  1 6QAM  in  System  View 

Figures  14  and  15  show  the  simulation  result  for  16QAM.  Again,  there  is 
excellent  agreement  between  the  simulation  results  and  the  analytical  results  in  equation 
(4.1). 


Performance  of  16QAM  in  AWGN 
with  no  channel  fading. 


Final  Value  Sink  8:  BER 
Loop  N  Time  Value 

10728  1.07350000000000  1. 3880788701 7429e-1 
20470  2. 04820000000000  9.52925481 227368e-2 
52G2G  5. 2G440000000000  5.734821 838201 56e-2 
107GG0  1 0.7G84000000000  2.72558781 843951  e-2 
27G925  27.  G955000000000  8.861 84385431 129e-3 
1104119  110.41 5500000000  1 . 81 3358381 23342e-3 
2446290  244.633200000000  1. 37091 324428854e-4 
3788461  378.850827000000  2.2351 8463742697e-6 


Figure  14.  The  simulation  of  16QAM  in  SystemView. 


Performance  of  simulated  16QAM  in  AWGN  without  coding 


Figure  15.  The  differences  of  performance  between  simulated  and  theoretical  16QAM. 
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c.  Simulation  of  1 6PSK  in  System  View 

Figures  16  and  17  show  the  simulation  result  for  16PSK  which  again 
compare  very  well  with  analytical  results  in  equation  (4.1). 


Performance  of  16PSK  in  AWGN 
with  no  channel  fading. 


Final  Value  Sink  18: 

Loop  N  Time  Value 

5807  591 .2000000000Q0e  3  1.G931933G2G8202e-1 
13528  1.35380000000000  1.3121G375803700e-1 
23695  2.37100000000000  9.83767830791 933e-2 
3838G  3.84060000000000  G. 807351 94009530e-2 
63947  6.39720000000000  3.9123630G729265e-2 
111362  11.1392000000000  2.1 09081 70582528e-2 
254397  25.4432000000000  6.991 34471 524253e-3 
925479  92.5520000000000  1.35750921 G45939e-3 
1596561  159.660863000000  1.31131515969732e-4 
2267643  226.769726000000  2.980261 72658483e-6 


Figure  16.  The  simulation  of  16PSK  in  SystemView. 


Performance  of  Simulated  16PSK  in  AWGN  without  coding 


Figure  17.  The  differences  of  performance  between  simulated  and  theoretical  16PSK. 
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B.  PERFORMANCE  OF  COHERENT  SYSTEMS  WITH  CODING 
1.  Convolutional  Codes  with  Hard  Decision  Decoding 

When  a  code  rate  r  =  k/  n  convolutional  code  is  employed,  the  probability  of  bit 
error  is  upper  bounded  by  [16] 

1  oo 

?,<7E«  (4.2) 

k  d-d 

a~a free 


where  dfree  is  the  free  distance  of  the  convolutional  code,  Bd  is  the  total  number  of 

information  bit  ones  on  all  weight  d  code  sequences,  where  d  represents  the  number  of 
code  bit  one  in  the  sequence,  and  Pd  is  the  probability  of  selecting  a  weight  d  output 
sequence  as  the  transmitted  code  sequence  assuming  the  all  zeros  sequence  was 
transmitted.  The  free  distance  of  the  convolutional  code  dfree  is  defined  as  the  weight  of 

the  minimum  weight  of  non-zero  code  sequence.  The  quantities  Bd  and  dfree  are 


parameters  of  the  convolutional  code  chosen,  and  Pd  is  determined  by  the  type  of 

modulation,  the  channel,  and  whether  hard  decision  decoding  (HDD)  or  soft  decision 
decoding  (SDD)  is  used.  For  HDD  [16] 


pt=i 


I 

i=(d+ 1)/2 

2 


f 


\l  j 


d-i 


(d  "1 

d 

Kdt  2) 

[p(i-j’)r"2+  Z 

'  i=l+d/2 

j 

P'(\-P) 


d-i 


for  d  odd 


for  d  even 


(4.3) 


where  p  is  the  probability  of  channel  bit  error. 

The  free  distances  and  the  first  eight  Bds  of  commonly  used  r  =  1  /  2 

convolutional  codes  are  listed  in  Table  2  [9],  and  the  corresponding  generator 
polynomials  are  listed  in  Table  3  [17].  The  generator  polynomials  are  a  concise  way  of 
specifying  the  connections  from  convolutional  encoder  shift  register  input  and  output  to 
the  coded  outputs.  For  example,  the  generator  polynomials  (131,171)  (expressed  in  octal) 
define  the  r  =  1/2  convolutional  encoder  shown  in  Figure  18.  [18] 
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Table  2.  Best  (maximum  free  distance)  rate  1/2,  convolutional  code  information  weight 
_ structure  (From  Ref.  9) _ 


V 

d  free 

B, 

O'  free 

B,  +, 

d  free  +1 

Bi  ? 

d  free 

B, 

“/ree+3 

B<tfree+4 

5,  +, 

dfree+5 

BJfree+6 

B,  7 

d  free  ' 

3 

5 

1 

4 

12 

32 

80 

192 

448 

1024 

4 

6 

2 

7 

18 

49 

130 

333 

836 

2069 

5 

7 

4 

12 

20 

72 

225 

500 

1324 

3680 

6 

8 

2 

36 

32 

62 

332 

701 

2342 

5503 

7 

10 

36 

0 

211 

0 

1404 

0 

11633 

0 

8 

10 

2 

22 

60 

148 

340 

1008 

2642 

6748 

9 

12 

33 

0 

281 

0 

2179 

0 

15035 

0 

Table  3.  Best  (maximum  free  distance)  rate  1/2,  convolutional  code  generators  (in  octal) 
_ (From  Ref.  17). _ 


V 

Generators 

dfree 

Bd 

u  free 

Bdfree+5 

Bdfree+6 

Bdfree+i 

3 

7,5 

5 

1 

192 

448 

1024 

4 

17,15 

6 

2 

333 

836 

2069 

5 

35,23 

7 

4 

500 

1324 

3680 

6 

75,53 

8 

2 

701 

2342 

5503 

7 

171,133 

10 

36 

0 

11633 

0 

8 

371,247 

10 

2 

1008 

2642 

6748 

9 

753,561 

12 

33 

0 

15035 

0 
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Figure  18.  The  r  =  1/2 ,  constraint  length  7  convolutional  encoder  defined  by  the  generator 

polynomials  (131,171).  (From  Ref.  18) 

2.  Simulation  of  Convolutional  Coded  BPSK  with  Hard  Decision 
Decoding 

Figures  19  and  20  show  the  simulation  results  for  BPSK  in  AWGN  with  r  =  1  /  2 , 
v  =  3  convolutional  source  coding  and  HDD,  which  compare  well  with  the  analytical 
upper  bound  in  equation  (4.2)  for  higher  SNR. 


Baseband  BPSK  in  AWGN  with  r  -  1  /  2  , 
v  -  3  convolutional  source  coding  and 
Hard  Decision  Decoding  (HDD). 


Final  Value  Sink  20:  coded  BER 

Loop  N  Time  Value 

1  231 3  2. 34G00000000000e-3  2. 1 G1  G947G8G98GGe-1 

2  GG32  G.  G98QQDQ0Q0Q0Q0e  3  1 . 1 57G75387821 25e  1 

3  1 4G55  1 4. 7540000000000e-3  G.  2320827G205908e-2 

4  289G3  29. 0950000000000e-3  3.49454850433324e-2 

5  71811  71.97G0000000000e-3  1.1GG915G08GG318e-2 

G  223730  223. 928000000000e-3  3. 291 2275G205G1  Oe-3 

7  1 003754  1 . 00398500000000  G.  41 00591 77GGG33e-4 

8  2G81 442  2.  G81 70G00000000  9. 9541 75031 352G7e-5 

9  43591 30  4. 35942750000000  1.19211 G770221 88e-5 


Figure  19.  The  simulation  of  baseband  BPSK  in  AWGN  with  r  =  1/2  ,  v  =  3  convolutional 

source  coding  and  HDD. 
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Performance  of  BPSK  in  AWGN  with  r  =  1/2,  v  =  3  convolutional  source  coding  and  Hard  Decision  Decoding 


Figure  20.  Performance  for  BPSK  with  r  =  1/2  ,  v  =  3  convolutional  source  coding  and 

HDD. 

3.  Simulation  of  Convolutional  Coded  BPSK  with  Hard  Decision 
Decoding  and  Barrage  Noise  Jamming 

The  primary  motivation  for  simulating  communication  systems  is  to  provide  an 
evaluation  tool  for  systems  that  are  difficult  to  analyze.  In  this  subsection,  this  principle  is 
illustrated  by  using  the  example  of  BPSK  with  r  =  1/2  ,  v  =  3  convolutional  coding 
where  the  signal  is  affected  by  barrage  noise  as  well  as  AWGN.  Since  the  barrage  noise 
is  just  more  AWGN,  we  expect  the  BER  performance  for  cases  where  El/N,«El/N0 

to  be  approximately  the  sum  as  for  the  corresponding  Eb  /  N0  in  the  AWGN  only  case,  as 
illustrated  in  Figure  20.  The  simulation  results  for  Eb  /  N0  =  13.35  dB  and  HDD  are 
shown  in  Figures  21  and  22  and  correspond  well  with  what  is  expected  when 

EbIN,«EbIN„. 
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Baseband  BPSK  (Eb/No=13.35  dB) 
in  AWGN  with  r  -  1  /  2  ,  v  =  3 
convolutional  source  coding, 
barrage  noise  jamming,  and 
Hard  Decision  Decoding  (HDD). 


Final  Value  Sink  14:  coded  BER 

Loop  N  Time  Value 

1  4804  4. 83700000000000e-3  2.081 598GG777G85e-1 

2  1 325G  1 3. 3220000000000e-3  1.1831 51 91 G70G1 1  e-1 

3  274G2  27. 5G1 0000000000e-3  7. 0392791 7781 21 9e-2 

4  5G45G  5G.  5880000000000e-3  3.44898944G09230e-2 

5  1 39333  1 39. 498000000000e-3  1 . 20660738202396e  2 
G  4G3808  4G4.  OOGOOOOOOOOOe  3  3. 081 901 53324G01  e-3 

7  1 983923  1 . 9841 5400000000  G.  578449G5G7G939e-4 

8  3GG1G11  3.  GG1 87500000000  9.1 792991 3070845e-5 

9  5339299  5. 33959G50000000  1 . 1 3251 0931 71 078e-5 

1 0  701 G987  7. 01 731 800000000  1.19211 G770221 88e-G 

11  8G94G75  8.G9503950000000  0.00000000000000e+0 


Figure  21.  The  simulation  of  baseband  BPSK  in  AWGN  with  r  =  l/2  ,  v  =  3  convolutional 

source  coding,  HDD,  and  barrage  noise  jamming. 


BPSK  performance  with  FEC,  HDD,  and  barrage  noise  jamming 


Figure  22.  Performance  for  BPSK  with  r  =  l/2  ,  v  =  3  convolutional  source  coding,  HDD, 

and  barrage  noise  jamming. 


28 


4.  Convolutional  Codes  with  Soft  Decision  Decoding 

For  soft  decision  detection  with  BPSK,  we  assume  that  the  receiver  is  a  maximal 
ratio  combiner. [7]  The  probability  of  bit  error  for  BPSK  in  AWGN  with  rate  r  =  1/2 
convolutional  source  coding  and  soft  decision  decoding  is  of  the  form  [15] 


1 

2(1  +  f) 

d 

pvn 

-brdCyh 

I'Ikvc 

2(1  +  ()  +  brh_ 

CAp 

2(l+()  +  brn_ 

(4.4) 


where  b  is  given  in  Table  1,  c  =  1.0  +  0.1^  is  empirically  obtained,  £  is  the  Ricean  factor, 
and  yh  is  the  average  Ebl  N0. 

5.  Simulation  of  Convolutional  Coded  BPSK  with  Soft  Decision 
Decoding 

Figures  23  and  24  show  the  simulation  results  for  BPSK  in  AWGN  with  r  =  1/2  , 
v  =  3  convolutional  source  coding  and  SDD,  which  compare  well  with  the  analytical 
upper  bound  in  equation  (4.2)  for  high  SNR.  For  Eb  / N0>  6  dB,  simulation  results  were 
not  obtained  due  to  the  excessive  runtime  required. 


RjR— 


kh  :  !  m  :  pm  : 

EB  laHm  i 


iBaseband  BPSK  in  AWGN  with  r-  1/2 
=  3  convolutional  source  coding  and 
II  Soft  Decision  Decoding  (SDD). 


Final  Value  Sink  20:  coded  BER 
Loop  N  Time 

1  3341  9.3740D0000QQ000e  3 

2  32519  32. 5850000000000e-3 

3  35520  95  61 90000000000e  3 

4  375716  375.  848000000000e- 

5  1 300536  1 . 300701 00000000 

6  3506224  3.58642200000000 

7  526391 2  5. 2641 4350000000 


Value 

2.35049056603774e-1 
8.71 0801 39372022e-2 
3.5756853396901 1e-2 
3  1.256281 4070351 8e-2 
4 17362270450751 e-3 
4.86381 3223571 98e-4 
6. 55B64223622032e  5 


Figure  23.  The  simulation  of  baseband  BPSK  in  AWGN  with  r  =  1/2  ,  v  =  3  convolutional 

source  coding  and  SDD. 
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Performance  of  BPSK  in  AWGN  with  r  =  1/2,  v  =  3  convolutional  source  coding  and  Soft  Decision  Decoding 


Figure  24.  Performance  for  BPSK  with  r  =  1/2  ,  v  =  3  convolutional  source  coding  and 

SDD. 

6.  Simulation  of  Convolutional  Coded  BPSK  with  Soft  Decision 
Decoding,  and  Barrage  Noise  Jamming 

Analogous  to  section  IV.B.3,  in  this  subsection  the  effect  of  both  barrage  noise 

and  AWGN  on  a  system  utilizing  BPSK  with  r  =  l/2  ,  v  =  3  convolutional  coding  is 

simulated,  only  here  SDD  is  used.  As  with  HDD,  simulation  results  are  obtained  for 

Eb/  N0=  13.35  dB.  We  except  the  BER  performance  for  cases  where  Ebl  NI  «  Eb  /  N0 

to  be  approximately  the  sum  as  for  the  corresponding  Eb  /  N0  in  the  AWGN  only  case,  as 
illustrated  in  Figure  24.  The  results  are  shown  in  Figures  25  and  26  and,  as  with  HDD, 
correspond  well  with  what  is  expected  when  Eb  /  N0 »  Eh  /  N : . 
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Baseband  BPSK  (Eb/No=13.35  dB) 
in  AWGN  with  r  -  1/2,  v  =  3 
convolutional  source  coding,  barrage 
noise  jamming,  and  Soft  Decision 
Decoding  (SDD). 


Final  Value  Sink  20:  coded  BER 

Loop  N  Time  Value 

1  1 0020  1 0. 0530000000000e  3  9. 9800399201 59G8e-2 

2  34G98  34. 7G40000000000e-3  4.0521 92235999G8e-2 

3  101941  1 02  040000000000e-3  1.487143G4320450e-2 

4  357211  357. 343000000000e-3  3.91 74207701 G492e-3 

5  1874024  1.87418900000000  G.59277049972541e-4 

G  3551712  3.55191000000000  9.1 792991 3070845e-5 

7  5229400  5. 229G31 50000000  7. 1 52700G21 331 2Ge-G 

8  6907088  G. 90735300000000  0.00000000000000e+0 

9  858477G  8.58507450000000  0.00000000000000e+0 


Figure  25.  The  simulation  of  baseband  BPSK  in  AWGN  with  r  =  1/2  ,  v  =  3  convolutional 

source  coding,  SDD,  and  barrage  noise  jamming. 


BPSK  performance  with  FEC,  SDD,  and  barrage  noise  jamming 


Figure  26.  Performance  for  BPSK  with  r  =  1/2  ,  v  =  3  convolutional  source  coding,  SDD, 

and  barrage  noise  jamming. 
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C.  FADING  CHANNELS 

1.  Performance  of  Coherent  System  with  Flat,  Slow  Fading  Channels 

For  BPSK  received  through  a  flat,  slow,  Rayleigh  fading  channel,  the  bit  error 
ratio  (BER)  is,  [15] 


and  yh  is  the  average  Ebl  N0. 

For  MPSK  and  Rayleigh  fading,  the  BER  is,  [15] 

qyb  sin 2{tt  I M) 
l  +  qyb  sin 2(7r/ M) 

For  MQAM  with  q  even  and  Rayleigh  fading,  the  BER  is,  [15] 

2(2'-l)  +  3  qrt/ 

For  MQAM  with  q  odd  (q  >  3)  and  Rayleigh  fading,  the  BER  is,  [15] 

3  qrb  _ 

2(2*-l)  +  3 qyb  t 

For  BPSK  and  Ricean  fading,  the  BER  is,  [15] 

exp 

where  c  =  1 .2  +  0.  \C,  is  empirically  obtained. 


n+C+ 1 


(4.6) 


(4.7) 


(4.8) 


(4.9) 
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For  MPSK  and  Ricean  fading,  the  BER  is,  [15] 

Pt =-_ - - —  expf  Z^EiM)£n_\  (4,10) 

q4nc[ \  +  £  +  q sin2 (n / M)/b]  ^1  +  ^  +  q sin2 (n / M)yb  J 

where  c  is  empirically  obtained. 

For  MQAM  with  q  even  and  Ricean  fading,  the  BER  is,  [15] 

p..  4(i-i/y'ixy-i)(4-+i >  f  ) 

‘  ?^[39n+2(2«-l)(f  +  l)]  {3qrb+2(2'-l)((  +  »J 

where  c  =  1.3  +  0.05^  is  empirically  obtained. 

For  MQAM  with  q  odd  and  Ricean  fading,  the  BER  is,  [15] 

_  ~3q^r>> - 1  (4.12) 

3qyb+2(V-\)^  +  \)) 

where  c  =  1.3  +  0.05^  is  empirically  obtained. 

2.  Simulation  of  Coherent  Systems  in  SystemView  with  Fading  Channels 
a.  Rayleigh  Fading 

(1)  BPSK.  The  performance  of  baseband  BPSK  transmitted  over  a 
slow,  flat  Rayleigh  fading  channel  is  shown  in  Figures  27  and  28.  The  results  of  the 
analytical  expression  in  equation  (4.5)  are  also  shown  in  Figure  28  and  shows  excellent 
agreement  with  the  simulated  results. 


4(l-l/2g/2)(^  +  1) 

q^[3qyb+2(2q-l)(C  +  l)] 
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Final  Value  Sink  10: 

Loop  N  Time 

1  235G2  23.5G20000000000e-3 

2  84792  84  /920000000000c  3 

3  284343  284.343000000000e- 

4  89202G  892.02G000000000e- 

5  2929494  2.92949400000000 

G  5890G12  5  89061200000000 

7  963391 0  9.63391 000000000 


Value 

G.3GG1 8283G771 07e-2 
2.331 3G4G2542742e-2 
3  7.51  G8753852398Ge-3 
3  2.4G839223740009e-3 
7.G03803321 0371 4e-4 
2.474G41 1 1 51 7227e-4 
8.1 062394201  G54Ge-5 


Figure  27.  Baseband  BPSK  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel. 


Performance  of  baseband  BPSK  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel 


Figure  28.  The  performance  for  simulated  and  approximate  baseband  BPSK  in  AWGN  over 

a  slow,  flat  Rayleigh  fading  channel. 


34 


(2)  16QAM.  Figures  29  and  30  show  the  simulation  results  for 
16QAM  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel.  There  is  excellent 
agreement  between  the  simulation  results  and  the  analytical  results  in  equation  (4.7). 


No  Coding  16QAM  in  AWGN 
over  a  slow,  flat  Rayleigh 
fading  channel 


Final  Value  Sink  8:  BER 
Loop  N  Time 

1  81 1 G  81 2. 200000000000e 

2  18548  1. 8581 0000000008 

3  32502  3.25200000000000 

4  45451  4.54750000000000 

5  58277  5.83070000000000 

6  71118  7.11540000000000 

7  83844  8.38880000000000 

8  87281  8.73280000000000 


Value 

3  9.451 785841 20983e-2 
4.1 4421 881 475342e-2 
1.40G6G760444507e-2 

4  583G51 G424751 7e-3 
I  52788388082506c  3 
4.771 44765721 820e-4 
1. 430956355831 15e-4 
4. 76985451 94371 6e-5 


Figure  29.  16QAM  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel. 


Performance  of  16QAM  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel 


Figure  30.  The  performance  for  simulated  and  approximate  16QAM  in  AWGN  over  a  slow, 

flat  Rayleigh  fading  channel. 
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(3)  16PSK.  Figures  31  and  32  show  the  simulation  results  for  16PSK 
in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel.  There  is  good  agreement  between 
the  simulation  results  and  the  analytical  results  in  equation  (4.6). 
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No  Coding  16PSK  in  AWGN 
over  a  slow,  flat  Rayleigh 
fading  channel 


Final  Value  Sink  18: 

Loop  N  Time  Value 

1  721  72.G000000000000e-3  1.07642G2G48008Ge-1 

2  1184  1 20. 400000000000e  3  5.5035773252G142e-2 

3  5873  588.8G0GGGGGGG80e-3  2.28258388495777e-2 

4  14718  1.47380000000000  7.729745204G9510e-3 

5  45434  4.54580000000000  2.4811 52781 75398e-3 

G  254233  25.4323000000000  8.588G057829945Ge-4 

7  588831  58. 88G7000000000  2.8G28G859433152e-4 

8  9253G9  92.5411310000000  9. 54289531 443840e-5 


Figure  31.  1 6PSK  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel. 


Performance  of  16PSK  in  AWGN  over  a  slow,  flat  Rayleigh  fading  channel 


Figure  32.  The  performance  for  simulated  and  approximate  16PSK  in  AWGN  over  a  slow, 

flat  Rayleigh  fading  channel. 
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b.  Ricean  Fading,  C,  -  4 

(1)  BPSK.  Figures  33  and  34  show  the  simulation  results  for 
baseband  BPSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  4 .  There  is 
good  agreement  between  the  simulation  results  and  the  analytical  results  in  equation  (4.9). 
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Final  Value  Sink  10: 

Loop  N  Time  Value 

1  32153  32.1 530000000000e-3  3.110129G9240817e-2 

2  255G80  255  680000000000c  3  5. 55308751  GG5926e-3 

3  15G2441  1.56244100000000  1. 02998397802701 e-3 

4  4917884  4.91788400000000  2.42292895453745e-4 

5  8273327  8.27332700000000  7.1 525641 94263G4e-5 

G  11G28770  1 1 .6287700000000  2.1457G925827909e-5 


Figure  33.  Baseband  BPSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  4 . 
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Performance  of  baseband  BPSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  4 


Figure  34.  The  performance  for  simulated  and  approximate  Baseband  BPSK  in  AWGN  over 

a  slow,  flat  Ricean  fading  channel  with  C,  =  4 . 


(2)  16QAM.  Figures  35  and  36  show  the  simulation  results  for 
16QAM  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  4  .  There  is  excellent 
agreement  between  the  simulation  results  and  the  analytical  results  in  equation  (4.1 1). 


Final  Value  Sink  8:  BER 

Loop  N  Time  Value 

1  11514  1.15200000000000  6.471 49460708783e-2 

2  32810  3.28220000000000  1.G9491525423729e-2 

3  58481  5.95090000000000  3.0G7484GG257GG9e-3 

4  89357  8.83810000000000  5.7257371 886G304e-4 

5  118723  11.9753000000000  1.43095G35583115e-4 

6  148013  14.8048000000000  4.76805834842881  e-5 


Figure  35.  16QAM  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  ^  =  4 
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Performance  of  baseband  16QAM  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  4 


Figure  36.  The  performance  for  simulated  and  approximate  16QAM  in  AWGN  over  a  slow, 

flat  Ricean  fading  channel  with  C,  =  4 . 

(3)  16PSK.  Figures  37  and  38  show  the  simulation  results  for  16PSK 
in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  C,  =  4  .  There  is  very  good 
agreement  between  the  simulation  results  and  the  analytical  results  in  equation  (4.10). 


No  Coding  16PSK  in  AWGN 
over  a  slow,  flat  Ricean 
fading  channel,  £  =  4. 
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Final  Value  Sink  18: 

Loop  N  Time  Value 

1  10171  1.01760000000000  7.9481 255861 8442e-2 

2  42470  4.24800000000000  3.698224852071 01  e-2 

3  210236  21.0252000000000  8.795411 0888661 6e-3 

4  378002  37.8024150000000  1.529636711 281 07e-3 

5  545768  54.5786300000000  2. 820734833041 26e-4 

6  713534  71.3568450000000  7.74888833255844e-5 


Figure  37.  16PSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  C,  -  4  . 
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Performance  of  baseband  16PSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  1,-4 


Figure  38.  The  performance  for  simulated  and  approximate  16PSK  in  AWGN  over  a  slow, 

flat  Ricean  fading  channel  with  £  =  4  . 


c.  Ricean  Fading,  C  =  \0 

(1)  BPSK.  Figures  39  and  40  display  the  simulation  results  for 
baseband  BPSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  C,  =  10.  There  is 
good  agreement  between  the  simulation  results  and  the  analytical  results. 
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Figure  39.  Baseband  BPSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  10 . 


Performance  of  baseband  BPSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  10 
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Figure  40.  The  performance  for  simulated  and  approximate  Baseband  BPSK  in  AWGN  over 

a  slow,  flat  Ricean  fading  channel  with  C,  =  10 . 


41 


(2)  16QAM.  Figures  41  and  42  display  the  simulation  results  for 
16QAM  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  <^  =  10.  There  is  good 


agreement  between  the  simulation  results  and  the  analytical  results. 


Figure  41 .  16QAM  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  10 . 


Performance  of  16QAM  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  ?  =  10 
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Figure  42.  The  performance  for  simulated  and  approximate  16QAM  in  AWGN  over  a  slow, 

flat  Ricean  fading  channel  with  C,  =  10 . 
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(3)  16PSK.  Figures  43  and  44  display  the  simulation  results  for 
16PSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  C,  =  10 .  There  is  a  good 


agreement  between  the  simulation  results  and  the  analytical  results. 


Final  Value  Sink  18: 

Loop  N  Time  Value 

1  G07  61.2000000000000e-3  6.471 49460708783e-2 

2  2462  247.200000000000e-3  2. 696871 62891 046e-2 

3  31412  3.1 4278888688808  3.14915545376467e-3 

4  1373583  137.360400000000  1. 430956355831 15e-4 

5  2715754  271.578127000000  5.503577325261 42e-6 

6  4057925  405.795854000000  5.9G057674540589e-7 


Figure  43.  16PSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  C,  =  10 . 


Performance  of  baseband  16PSK  in  AWGN  over  a  slow,  flat  Ricean  fading  channel  with  £  =  10 


Figure  44.  The  performance  for  simulated  and  approximate  16PSK  in  AWGN  over  a  slow, 

flat  Ricean  fading  channel  with  <^  =  10. 
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C.  CHAPTER  SUMMARY 

In  this  chapter,  the  bit  error  ratio  performance  of  various  types  of  modulated 
signals,  with  and  without  error  correction  coding  and  with  and  without  channel  fading, 
were  obtained  by  simulation  and  the  results  compared  with  those  obtained  analytically 
when  possible.  In  the  next  chapter,  conclusions  are  drawn  and  suggestions  for  future 
work  made. 
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V.  CONCLUSION  AND  FUTURE  WORK 


A.  CONCLUSION 

Using  SystemView,  we  simulated  the  effect  of  various  types  of  fading  channels. 
Two  major  problems  were  found  during  the  modeling  and  simulation.  First,  the  total 
number  of  computations  is  restricted  on  each  simulation  loop,  which  limits  the  degree  of 
accuracy  for  high  SNR.  Next,  too  much  CPU  resource  is  required  for  each  simulation, 
and  it  is  time  consuming  to  accurately  simulate  performance  for  high  SNR.  Presumably, 
using  the  commercial  version  instead  of  the  student  version  would  alleviate  these 
difficulties. 

Except  for  these  two  problems,  generally  speaking,  SystemView  provides  a  good 
modeling  and  simulation  tool  for  communication  systems.  The  performance  of  coherent 
systems  with  no  coding,  the  performance  of  coherent  systems  with  coding,  and  the 
performance  of  coherent  systems  with  no  coding  and  fading  channels  all  agree  very  well 
with  analytical  results.  Most  importantly,  the  performance  of  communication  systems 
affected  by  narrowband  noise  and  other  channel  impairments  that  are  difficult  to  analyze 
are  easily  simulated  in  SystemView.  Since  the  simulations  developed  have  been  verified 
for  the  AWGN  case,  they  can  now  be  modified  to  examine  cases  that  are  not  possible  to 
evaluate  by  analysis. 


B.  FUTURE  WORK 

With  rapid  progress  in  hardware  and  software,  the  simulation  process  will  be 
significantly  improved.  The  utility  of  this  simulation  tool  for  evaluating  difficult  to 
analyze  systems  has  been  demonstrated,  and  SystemView  can  be  used  to  evaluate  any 
number  of  modem  communication  systems.  The  combination  of  error  correction  coding 
and  modulation  specified  by  the  various  digital  standards,  such  as  IEEE  802.1  lg,  should 
be  examined  for  fading  channels  and  interference  conditions  other  than  simple  AWGN  to 
foster  improved  understanding  of  how  these  popular  communications  wireless  standards 
perform  in  various  radio  environments.  For  example,  the  barrage  noise  jamming 
considered  in  Chapter  IV  is  simply  additional  AWGN.  It  would  be  valuable  to  re- 
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examine  this  case,  but  for  a  barrage  noise  jammer  implemented  by  frequency  modulating 
a  voltage  controlled  power  oscillator  using  bandlimited  AWGN  or  bandlimited  AWGN 
plus  a  periodic  waveform,  as  is  typical  of  real  military  electronic  warfare  systems  [19]. 
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